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Abstract 



o 

r\ ' We analyse the dynamical density fluctuation spectra for cuprates starting 

from the t~ J model in a slave-boson 1/N representation. The results obtained 
are consistent with diagonalization studies and show novel low-energy structure 
on the energy scale J + St due to the correlated motion of holes in a RVB 
spin liquid. The low-energy response implies an anomalous renormalization 
of several phonon modes. Here we discuss the renormalization of the highest 



p^ ' breathing phonons in La2-a;Sri:Cu04. 

O, 

1. INTRODUCTION 

C^ I High-temperature superconductors are doped Mott-Hubbard insulators, there- 

^ ■ fore the low-energy density response is proportional to the doping, due to tran- 

sitions in the lower Hubbard band. In one dimension the Hubbard physics is 
^ ■ characterized by charge and spin separation, which implies that the density 

P ■ response is spinless fermion like, i.e., showing vanishing excitation energy at 

Akp- Exact diagonalization studies [1,2] forthet-J model have revealed that the 

dynamical density response A^(q, w) for the 2D model relevant for the cuprates 

K> ', is very different from the ID case. On the other hand these calculations show 

^ [ several features also unexpected from the point of view of weakly correlated 

fermion systems: (i) a strong suppression of low energy 2kp scattering in the 
density response, (ii) a broad incoherent peak whose shape is rather insensitve 
to hole concentration and exchange interaction J, (iii) a very different form of 
A^(q, u) compared to the spin response function 5(q, a;), which share common 
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features in usual fermionic systems. Finite temperature diagonalization studies 
[3] show only weak temperature dependence for T < 0.3t even at low energy. 

While considerable analytical work has been done to explain the spin re- 
sponse of the t-J model only few authors analysed N{q,uj). Wang et al.[4] 
studied collective excitations in the density channel and found sharp peaks at 
large momenta corresponding to free bosons. Similar results were obtained by 
Gehlhoff and Zeyher[5] using the X-operator formalism. Lee et al[6] consid- 
ered a model of bosons in a fluctuating gauge field and found a broad incoherent 
density fluctuation spectrum at finite temperature, due to the coupling of bosons 
to a quasistatic disordered gauge field. 

Starting from a slave-boson representation we show that the essential features 
observed in the numerical studies can be obtained in the framework of the 
Fermi-liquid phase of the t-J model at zero temperature [7]. Our main findings 
are: (i) at low momenta the main effect of strong correlations is to transfer 
spectral weight from particle-hole excitations into a pronounced collective 
mode. Because of the strong damping of this mode (linear in q) due to the 
coupling to the spinon particle-hole continuum, this collective excitation is 
qualitatively different from a sound mode, (ii) At large momenta we find a 
strict similarity of iV(q, cj) with the spectral function of a single hole moving 
in a uniform RVB spinon background. In this regime iV(q, a;) consists of a 
broad peak at high energy whose origin is the fast, incoherent motion of bare 
holes, (iii) The polaronic nature of dressed holes leads to the formation of a 
second peak at lower energy, which is more pronounced in (vr, 0) direction in 
agreement with diagonalization studies[l, 2]. 

The anomalous renormalization of certain phonon modes as observed in 
inelastic neutron scattering provides a sensitive test of the pecularities of the 
low-energy density response. In this contribution we shall analyse the strong, 
doping-dependent renormalization of the highest breathing phonon modes 
which is a generic feature in the high-Tc compounds. 

2. SLAVE BOSON THEORY OF DENSITY RESPONSE 

Following Kotliar and Liu[8] and Wang et al [4] we start from the A^- 
component generalization ofthe slave-boson t-J Hamiltonian, Htj = Ht+Hj, 
which is obtained by replacing the constrained electron creation operators 

^* = -f E {fthp^f,. + h.C.\ (1.1) 

<i,j>(T 

HJ = ^ E ftJ^<^'fta'fU^-Kh^){l-h+h,), (1.2) 
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Figure 1.1 Comparison of N{q,ij) obtained by the slave-boson theory (solid lines) with 
diagonalization data for a periodic 4x4 cluster with J/t = 0.4 and doping 5 = 0.25. The 
dashed line in the (tt, tt) spectrum indicates the (5-function collective peak obtained when polaron 
effects are neglected. Right figure shows the weak J-dependence for the (tt, tt) spectrum. 



where fj^^ is a fermionic (spinon) operator, a = 1, • • • , A^ is the fermionic 
flavor index, and hi denotes the bosonic holes. These operators obey standard 
commutation rules, yet the number of these auxiliary particles must obey the 
constraint J2a fi^ficr + hfhi = N/2. The original t-J model is recovered for 
N = 2. 

The slave boson parametrization provides a straightforward description of 
the strong suppression of density fluctuations of constrained electrons through 
the representation of the density response in terms of a dilute gas of bosons. 
A common treatment of model (1) is the density-phase representation ("ra- 
dial" gauge[9]) of the bosonic operator hi = Tj exp(i^j) with the subsequent 
1/A^-expansion around the Fermi-liquid saddle point. While this gauge is par- 
ticularly useful to study the low energy and momentum properties, it is not 
very convenient for the study of the density response in the full uj and q space. 
Formally the latter follows in the radial gauge from the fluctuations of rf. 
If one considers for example convolution type bubble diagrams, one realizes 
that their contribution to the static structure factor is correctly of order 1/N, 



but is not proportional to the density of holes 6 as it should be. According 
to Arrigoni et al[10] such unphysical results originate from a large negative 
pole contribution in the (r_qrq)i^ Green's function of the real field r, which 
is hard to control by a perturbative treatment of phase fluctuations. We follow 
therefore Popov[ll] using the density-phase treatment only for small momenta 
q < qo, while keeping the original particle-hole representation of the density 
operator, 6+6, at large momenta. More precisely hi = ri exp(i^j) + 6j, where 
bi = I]|q|>g,) ^q exp(iqRj). The cutoff qo is introduced dividing "slow" (col- 
lective) variables represented by r and 9 from "fast" (single-particle) degrees 
of freedom. As explained by Popov[ll] this "mixed" gauge is particularly 
useful for finite temperature studies to control infrared divergences. We start 
formally with "mixed" gauge and keep only terms of order 6 and 1/A^ in the 
bosonic self energies. In this approximation our zero temperature calculations 
become quite straightforward: The cutoff qo < 6 actually does not enter in the 
results and we arrive finally at the Bogoliubov theory for a dilute gas of bosons 
moving in a fluctuating spinon background. 




Figure 1.2 Density fluctuation spectra N{q, a;) for J/t = 0.4 and 5 = 0.15 along (n, 0) (top) 
and (n, n) directions (bottom). Energy in units of t. 
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The Lagrangian corresponding to the model (1) is then given by (the sum- 
mation over a is impUed) 

i 

+^T.>^^{ft.h<r + {n + bt){n + b,)-^), (1.3) 



2t 
^* = -]7 E ftafjAbpi + TiVj + rjbi + b+n) + h.c. (1.4) 

Here the A field is introduced to enforce the constraint, and jif, m, are fixed 
by the particle number equations (nj) = ^(1 ~ ^) ^rid {rf + bfbi) = y^' 
respectively. The uniform mean field solution r^ = ro \/N/2 leads in the large 
N limit to the renormalized narrow fermionic spectrum ^^ = —zi^k — ^/» 
with i = Jx + tS, 7k = ^(cos ka: + cos ky), x = T,a{fi^fj^)/N' and z = 4 
the number of nearest neighbors. In the N = oo limit Xoo — 2/7r^ is given 
by that of free fermions, while for the original t-J model its value should 
be larger[12] due to Gutzwiller projection. In the following x = fxoo will 
be used. Distinct from the finite-temperature gauge-field theory of Nagaosa 
and Lee[13] the bond-order phase fluctuations acquire a characteristic energy 
scale in this approach[4], and the fermionic ("spinon") excitations can be 
identified with Fermi-liquid quasiparticles. The mean field spectrum of bosons 
is ujq = 2zxt{^ — 7q)- Thus the effective mass of holes m^ oc 1/t is much 
smaller than that of the spinous. 

Due to the diluteness of the bosonic subsystem, 5 <^1, the density correla- 
tion function Xq,tj = {Sn^6n^)ciui is mainly given by the condensate induced 
part which is represented by the Green's function ((6+ + 6_q)(6q + 6lq))a; for 
q > qo, and 2(r_qrq)(^ for q < qo, respectively: 

Xqc. = -J^o{{(K + ^-q)(^q + b-ci))q>qo + 2(r_qrq)g<,„) . (1.5) 

The 1/A^ self-energy corrections to these functions are calculated in a conven- 
tional way [9, 8] expanding rj = (ro\/iV+(5r)j)/\/2 and considering Gaussian 
fluctuations around the mean field solution. Neglecting all terms of order 6/N 
and q^/N , only one relevant 1/A^ contribution remains which corresponds to 
the dressing of the slave-boson Green's function by spinon particle-hole excita- 
tions. Within this approximation and at zero temperature no divergences occur 
at low momenta, thus one can take the limit qo -^ 0. The final result for the 
dynamic stucture factor (normalized by the hole density) is: 

iVq,c. = ^/m((L^qa + 41/^) - ^ib)/D^,^) , (1.6) 



The origin of the contribution 

4i = *S(^Hip2 - "3),,., (1.8) 

n_^_^^,^ nteJ_n&^^ ._ (1.9) 

is the indirect interaction of bosons via the spinon band due to the hopping term 
(which gives U^ in (4)) and due to the coupling to spinons via the constraint 
field A. The latter channel provides a repulsion between bosons, making 
S^^' {uj = 0) positive and therefore ensuring the stability of the uniform mean- 
field solution. The 1/A^ self energies S^^'^' and A^^'^' are essentially a 
single boson property. They are given by the symmetric and antisymmetric 
combinations (with respect to u; + iO+ -^ —to — iO^) of the self energy 

^S/f'^ = 4 E (^t7k'-q)=^G'0+k-k'('^ + Ck - Ck')- (1-10) 

|k|<fcF<|k'| 

Here G'^{uj) = {uj — lo^ — Sq,^ + fib)^^ is the Green's function for a single 
slave boson moving in a uniform RVB background. Although in the context of 
1/A^ theory the G^ function in (5) should be considered as a free propagator, 
we shall use here the selfconsistent polaron picture for a single hole[14]. This 
is crucial when comparing the theory for A^ = 2 with diagonalization studies. 
Finally, the constants a and //^ in (3) are given by (1 — trQ/t) and S^'^^^^ {uj = 
q = 0), respectively. The parameter Tq in Eq. 1.8, which formally corresponds 
to the condensate fraction in our theory, is determined selfconsistently from 
rQ = 5 — Z]q3^o ^q- "^^^ momentum distribution n^ = (b^bq) is calculated 
from the corresponding bosonic Green's function for finite hole-density. 

In the small lo, q limit A^(q, to) (3) is mainly controlled by the interaction of 
bosons represented by the S^^' term (oc rg), while the internal polaron structure 
of the boson determined by S^^'^' is less important. A^(q, ut) consistsof a weak 
spinon particle-hole continuum with cutoff oc vpq, and a very pronounced linear 
collective mode which nearly exhausts the sum rule. The velocity of this mode 
is always somewhat smaller than the spinon Fermi velocity, Vs < vp — zi, 
which implies a strong damping oc to (or q) of this mode (Fig. 1.2). 

The density response A^(q, w) at large momenta, q > 6, which we can 
compare with diagonalization results, is dominated by the properties of a single 
boson selfenergy S^^'^\ The calculated density response of the t-J model has 
three characteristic features on different energy scales: (i) The main spectral 
weight of the excitations at large momenta is located in an energy region of 
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order of several t. This high energy peak is very broad and incoherent as a result 
of the strong coupling of bosons to low-energy spin excitations. The position of 
this peak and its shape are rather insensitive to the ratio J/t < 1 in agreement 
with conclusions of [1, 2]. This is simply due to the fact that the high-energy 
properties of the t-J model are controlled by t. (ii) The theory predicts also 
a second peak at lower energy (Fig. 1.2) which is more pronounced in the 
direction (vr, 0), while its weight is strongly suppressed for q near (tt, vr). The 
origin of this excitation is due to the formation of a polaron-like band of dressed 
bosons. The relative weight of this contribution increases with J as a result of 
the increasing spinon bandwidth, (iii) In addition there is the spinon particle- 
hole continuum which is generated by S^^^ with relatively small weight (oc 5"^). 
At (tTjO) the high energy cutoff of the spinon continuum is at z{xJ + 5t) 
(Fig. 1.2), while the polaron peak is at about {xJ + ^i)- 

We note that the polaron peak in A^(q, uj) is in the same energy range as the 
high energy phonons in cuprates. 

3. RENORMALIZATION OF BREATHING PHONONS 

Inelastic neutron scattering experiments on high- Tc superconductors have 
shown that in particular the highest energy longitudinal optic phonons near 
(vr, 0) soften and broaden strongly as holes are doped into the insulating parent 
compound. Whereas the corresponding breathing mode at (vr, vr) shows much 
smaller softening and no anomalous broadening. This effect seems to be 
generic for cuprates and detailed neutron scattering studies have been reported 
for La2-xSr^Cu04 [16, 17, 18] and YBa2Cu306+x [18]. 



Cu O 

a) '^ b) o 

Figure 1.3 Atomic displacements of oxygen ions (a) for tlie q = (tt, tt) breatliing plionon and 
(b) for the (tt, 0) lialf-breatliing mode. 



The renormalization of these phonons can be calculated in the framework 
of the t-J model, since these modes modulate the energy of a hole in a Zhang- 
Rice singlet state, and therefore couple directly to the density of doped holes. 

Expanding the Zhang-Rice energy Ezr = 8-^ with respect to the oxygen 
displacements u^, Vy (see Fig. 1.3) of the four 0-neighbors of the Cu-hole 



yields the linear electron-phonon coupling 



H,- 



e—ph 



9^{u'^ 



U_^ + Vy 



v%)hth 



y)"'i "'I- 



(1.11) 



—7/2 

We assume that the resonance integral obeys the Harrison relation tp^ oc Tq , 
where tq is the Cu-0 distance, and obtain g = lEzn/'irQ, i.e., g « 4eV/X. 
The lattice part of the Hamiltonian is determined by the force constant K k, 
25eV/X^ for the longitudinal 0-motion. Due to the structure of H^-ph the 
breathing modes couple directly to Xq,w- 

We have studied the renormalization of the phonon Green's functions along 
(vr, 0) and (vr, vr) directions 



jjph 



^qfi 



UJ^ 



<o(l 



Q^qXq,! 



(1.12) 



where Uq^ is the bare phonon frequency, i.e. measured in the undoped parent 

compound, and Qq = ^(sin^ qx/2 + sin^ '7j//2). Based on the parameters 
of the pd-model we estimate for the dimensionless coupling constant ^ = 

g'^/ztK ~ 0.3 - 0.5. 
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Figure 1.4 Calculated phonon spectral function Bph for (tt, 0) and (tt, tt) breathing phonons 
for 5 — 0.15, i = 0.4 eV, J = 0.12 eV and ^ = 0.25 (solid lines) compared to undoped system 
(dash-dotted lines) and the inelastic neutron scattering data for Lai.85Sro.i5Cu04. 



Figurel.4 shows the strong renormalization of the (vr, 0) half-breathing mode 
for Lai.85Sro.i5Cu04 with a twice as large shift as for the (vTjVt) breathing 
phonon. The large damping of the (tt, 0) phonon results from the hybridization 
with the large polaron peak in A'^(q, w) at this momentum[19] and is consistent 
with the experimental data[17]. The phonon energies of the undoped parent 
compound uoq^ = 80(90) meV for (vr,0) and {tt,!:), respectively, are taken 
from Ref.[16]. The strong doping dependence of this effect is shown in Fig. 1.5. 
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Figure 1.5 Doping dependence of low-energy density response at (n, 0) (solid liness). As a 
consequence of the scaling of the polaron structure ex (xJ + St) there is a strong change in the 
renormalization and damping of the (tt, 0) half-breathing phonon, which is at ujo = 0.2i in the 
undoped system. 
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Figure 1.6 Fano structure in fully renormalized N{q, uj) (solid line) due to the coupling to the 
(tt, 0) phonon. The bare A''(q, u>) and the phonon spectral function are indicated by dashed and 
dotted lines, respectively (parameters as in previous figure). 



Finally we have studied the changes of the density response Xq,tj due to the 
additional coupling to the breathing phonon modes. This effect is displayed in 
Fig. 1.6, which shows a rather strong Fano structure in A^(q, uj) for q = (vr, 0). 
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4. SUMMARY 

We have outlined a 1/N slave-boson theory for the density response of the 
t-J model, which explains the data obtained by exact diagonalization. We 
demonstrated that the predicted low energy polaron structure in the density re- 
sponse, which is particularly pronounced along (vr, 0), explains the anomalous 
doping induced line width and shift of the longitudinal planar (vr, 0) phonon. 
The energy of the polaron peak is determined by the spinon energy scale, there- 
fore we predict a nontrivial doping dependence for the phonon renormalization. 
In that respect further neutron scattering studies of the doping dependence of 
phonons would provide a sensitive test for the low energy density response as 
well as for the spin structure in the different doping regimes. 
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